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Abstract. Sometimes we obtain attractive results when associating 
facts to simple elements. The goal of this work is to introduce a possible 
alternative in the study of the dynamics of rational maps. In this study 
we use the family of maps f(x) = fyEfj making some associations with 
1 -a 
1 -b 



the matrix A — 



of its coefficients. 



1. Introduction 

The main goal of this article is to present an alternative tool to study 
the dynamics of a real rational function, using results from the Numerical 
Range Theory, and has two main parts. The first, comprising Sections 2-3, 
is concerned to the adequation of the Numerical Range Theory to the data 
provided from rational maps. As described in Milnor[8j, each map /, in the 
space Rat 2 , can be expressed as a ratio 

(Ml I I I t ~ 
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p(z) 



a z z + a\z + a 2 



q(z) b z 2 + b\z + b 2 ' 

where ao and bo are not both zero and p(z), q(z) have no common root. 
Milnor[8j states that we can obtain a roughly description of the topology 
of this space Rat 2 that can be identified with the Zariski open subset of 
complex projective 5-space consisting of all points 

(a : ai : a 2 : b : b\ : b 2 ) G CP 5 , 

for which the resultant 
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is non-zero. Taking z = x + iO and ao = 1, bo = 1, ai = —a+iO, b 2 
a\ = 0, b\ = 0, with x, a, b real numbers, we obtain 



-b+iO, 
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res(p,q) = det B. This matrix B is associated to the real rational map 

- 2 - a) / (t 2 
our results associated to the matrix B. 

i 



f(x) = (x 2 — a) I (x 2 — 6). The map / will be the one that we will use in 
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The second part of this article, comprising Sections 4-5, shows how can 
we apply the Numerical Range Theory to the dynamics of the map /, estab- 
lishing the relation between some partitions of an ellipse and the symbolic 
space generated by the partition of the domain of / in real intervals. More- 
over, we launch a conjecture that could be a path to generate, in the future, 
an extension of the usual symbolic space applied to rational maps, allowing 
us to describe much better the dynamics of this maps. 

2. Numerical Range Theory 

The classical numerical range of a square matrix M n , with complex num- 
bers elements, is the set W(M n ) = {u*M n u, u £ S{C n )} , with S(C n ) the 
unit sphere, u is a vector in C n and u* is the transpose conjugate of u. 
The numerical range W(M n ) also can be defined as the the image of the 
Rayleight quotient RM n (u) = u* M n u/u*u, The set W{M n ) is closed 

and limited, and it is also a subset of the Gaussian C plane. Toeplitz[10j and 
Hausdorff[2j proved that W(M n ) is a convex region. From Kippenhahn[5j, 
the boundary of the numerical range, dW(M n ) is a piecewise algebraic curve. 
In the particular case of a square matrix M2, with eigenvalues Ai, A2, W{M2) 
is a subset limited by an ellipse with foci in Ai and A2, result known as El- 
liptical Range Theorem, see Li[6|. 

Proposition 2.1. // H Mn = (M n + M*) /2 and S Mn = (M n - M*) /2 are 

the Hermitian and skew- Hermitian parts of M n , respectively, thenHe(W(M n )) = 
W(H Mn ) andlm(W(M n )) = W(S Mn ). 

Proof. The proof can be found in Melo[7j. □ 

Theorem 2.1. To every complex matrix M n = Hm„ + Sm„ through the 
equation kM n (ctx, 012, CK3) = &et{a\HM n — ioi^SMn + ct^I n ) = is associated 
a curve of class n in homogeneous line coordinates in the complex plane. 
The convex hull of this curve is the numerical range of the matrix M n . 

Proof. Adapting the proof in Kippenhahn[5] we have the desired result. □ 

3. Merging f(x) in W(M n ) 

Hwa-Long Gau[lj states that we can obtain from a 4 x 4 matrix an ellip- 
tical numerical range, thus we could use B as defined in section [TJ but we 
can simplify our results if we use a smaller matrix, A2, trough a result in 

linear algebra. The new matrix A2 = ^ | ^ j will produce equivalent 

results as the obtained from B. 

Lemma 3.1. The matrix B is unitary decomposable in 

( A 2 2 \ 
V ® 2 M J 

a block diagonal matrix. 

Proof. In order to prove this result it is sufficient to find an unitary matrix 
E such that 



E*BE = 
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With some computation we can see that 

/ 1 \ 

10 

10 

\ 1 J 

□ 

Proposition 3.1. W{A 2 ) = W{B) 

Proof. By lemma I3TT1 E*BE = A 2 ®A 2 and using the properties of numerical 
range we have 

W(E*BE) = convex hull{W(A 2 ) U W(A 2 )} = convex hull{W(A 2 )}. 

The convex hull of a convex set is itself, then W(E*BE) = W{A2). But 
the numerical range of B is invariant under unitary transformations, see 
Kippenhahn[5], so W(B) = W(A 2 ). □ 

In our study we use f(x) = (x 2 — a)/ (x 2 — b) , with a > 0, b > and a > b. 
Such map takes all real axis with exceptions ±y/b on (— oo,l) U [|,+oo). 
With A = M\{±^}x (-oo,l) U [f,+oo) we can define the graphic of /, 
graph(f), as the pair (x, f(x)) e A and 9 : R\{±y/b} — > A. 

Definition 3.1. Let C = {v G C 2 : v = (x, if(x)} and 

* = \ z £ C : z = — ,v / 

[ v*v 

We define V : A — ► C as (x, f(x)) i — > (x, if(x)) and E : C — > 

By definition 13.11 the image of (x, f(x)) is z = (v*A 2 v) /v*v by H o V. 

Proposition 3.2. ^ c W(A 2 ) 

Proof. By the definition of VP and VF(t42), using Rayleight quotient, the result 
follows. □ 

From proposition 13.21 we know that \P is a subset of VK(A2) C C, and 
using definition 13. 1\ we can calculate the elements z S f , and as they were 
defined, they will become a function of x. After some calculations we have 

-a 2 b + (2a + b) bx 2 - 3bx 4 + x 6 + i (a + 1) (-abx + (6 + a)x 3 - x 5 ) 
Z(X) = a 2 + (b 2 - 2a) x 2 + (1 - 26) x 4 + x 6 ' 

So, the z 6 f is a function such that z[x) = g{x) + ih(x), with Some 
elementary calculus show us that g(x) and h(x) are real rational continuous 
functions in R, therefore z(x) is continuous in C. We call some attention to 
the fact that z(±Vb) and z(oo) exists and are well defined in C. 



Observation 3.1. We have z(x + 1) = z(x) for 
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Let X2 and X12 be the values where f(x) = 0. If we calculate z(0,a/b); 
z(x2,0); z(xi2,0); z(x,x); z(x,—x) we obtain four different points of \&. 
With some elementary algebra we calculated the ellipse that contain this 
four points. This ellipse is 

n = < ± V- + 9 = 1, (x, y) e M 2 , a > b, a > 0, b > } , 

with 1 + 6 and 1 + a the minor and major axis length of O, respectively. 

Moreover, when we use all points (x, f(x)) they will fall in f2 under 
transformation by z. 

Lemma 3.2. Let z(x) = g(x) +ih{x), then the pair (Re(z), Im(z)) satisfies 
O. 

Proof. We obtain this result replacing in the equation of 0, x by Re(z) and 
y by Im(z). □ 

Proposition 3.3. If Si E \P £/zere are, a£ Zeasi one Xj such that z{xi) = Si. 

Proof. Since z{x) = g{x) + ih{x) is a continuous function in C and by the 
lemma the result follows. □ 



Then we conclude that ^ can be represented by the ellipse with equation 

n. 

Since is constructed in the space R 2 and this space is isomorphic to C, 
when we refer to an element z G it can understood has a vector in M? or 
a complex number in the plane C. 

There are relations between the functions / and g that we can observe, 
described in the following lemmas. The proofs are omitted because they 
result from straight calculus. 

Lemma 3.3. If xq is a zero of f{x), then g(xo) is a relative maximum of 

g(x). 

Lemma 3.4. If xo is a relative minimum of f(x) or xo is a discontinuity 
value of f{x), then g(xo) is a relative minimum of g(x). 

There are similar relations between h(x) and /(x). 
Follows some results relating /(x) to 17. 

Lemma 3.5. Let /(xo) = ±xo, then E o V(xq, /(xo)) is vertex of Q. 

Proof. If /(xo) = xo, by V we have (xq, ixq). So 

x -ix ) ( 



r.((x , tx )) = t ^ = — i- 

x 



ix 



And if we look at the equation of 0, we see that (^p, — ^2^) is a vertex 
of tt. 

If /(xo) = — xo we obtain another vertex of O in a similar way, which is 
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Lemma 3.6. The discontinuities of f(x) and the values where f(x) has a 
minimum are transformed by H o V o in the vertex (—6,0) of Q,, and the 
roots of f(x) and the oo are transformed by S o V o 9 in the vertex (1,0) of 

n. 

Proof. Since z(x) = g(x)+ih(x) = v* A2V /v*v, v = (x, if(x)) is a continuous 
function in C, we have z(±Vb) = —b and z(oo) = 1. □ 



4. Partitions of O 



Let x\, x<2, X5, xq, X7, xs, xq, X12, X13 be the solutions of g'(x) = 0. By 
lemma 13,41 and lemma 13.31 and considering the order of real axis, we will 
have X2 and xu as zeros of f{x); X5 and xg as the discontinuities of f(x) 
and x-j = 0. All this values have image from z(x), lemma [3T2| including 
the infinity, being related by z{x2) = z{x\2) = 21(00) and equal to vertex 
(1,0) in Q, see lemma [3761, z(x§ ) = z(xg) = z(0) and equal to vertex (—6,0) 
in Q. Related to the real axis, z(xi) is symmetric to z(xis) and z(xq) is 
symmetric to z(xs) in $7. Where are the missing X3, X4, xio, x\\ ? They 



will be the values such that z(x\) 
and z(x 13 ) = z(x 3 ). 

Using this special values Xi, i = 
define a partition function pa, as 



z(x n ), z(x 6 ) = z(x w ), z(x 8 ) = z(x 4 ) 
L, .., 13, with order Xi < a?i+i, we can 



pa{x 




if x < x\ 

if Xi-i < x < Xi with 2 < % < 13 
if x > X13 



Now we will create partitions in Q using the images z(xi), i = 1, .., 13 in 
£1. Here, we ask attention for one particular aspect of z, see proposition 13.31 
Some intervals Ii will be transformed in the same arc of Vl. The only thing 
that will distinguish them is the orientation and the origin of its end points. 

Definition 4.1. Let Si = z{xi), we define arc(Si,Si+i) as the arc of Q 

starting at Si and ending at S'i+i, with counterclockwise orientation. 



We define pan, a partition function, as 

pa Q (w) 



wee 



J\, if w G arc(z(oo), z{x\)) 
Ji, if w £ arc{z(xi-\), z{xi)) with 2 < i < 13 
J14, if w £ arc(z(xi4), z(oo)) 



The functions pa(x) and pan(w) are related by 



zili 



Ji if i = l 
-Ji if 2 < i < 6 

if i = 7 or i = 8 
-Ji if 9 < i < 13 
J14 »/ * = 14 



thus we can build a matrix T14 of the transformation pon(^(pa)), 



N 7 O r 
O7 iV 7 
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with 



N 7 



1 
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Il4, 


deb(T\4) = 


1, 


and 



matrix. 



5. Dynamics of f(x) 

Now we have a new tool to study the dynamics of / using a symbolic space. 
Using n, to study the behavior of / we will have the same advantages that 
we would have when studying the behavior of second degree polynomials 
functions in the unit circle. 

If we define a symbolic space using the partitions created by the function 
pa in the real axis we will have the problem of dealing with the discontinuities 
of the function / and the infinity itself. So, profiting that z is a continuous 
complex function in C U {oo} this problem will vanish. 

We can build two distinct symbolic spaces. The first will be the classi- 
cal association between the intervals produced by pa in the real axis, see 
MilnorJS] for further reference, using the domain of the function /, and con- 
sidering an alphabet A with designations Jj for each interval, we will have 
a symbolic space E c = A^. The second will be constructed as we consider 
the alphabet B = {J±, J14}, and the set E = £> N of symbolic sequences 
on the elements of B, introducing the map spa : R U {00} — > B. 

Conjecture 5.1. The symbolic dynamics of f does not change if we use £ 
instead of £ c . 

Both spaces are connected by the transformation matrix T14 and doing 
some calculus in matrix algebra, since this matrix is an involutary matrix, we 
could get the result. All computations in our work points in that direction. 
But we are still working in a suitable proof of this result. Moreover, £ will 
work as an extension of E c . 

It means that we can identify the periodic orbits in the same values of a 
and b as we use both spaces E and E c . For example for the values a = 4.01 
and b = 2.5 the critical orbit of / is periodic in both spaces, such as all the 
others values of periodicity found in our research. But they are many new 
sequences in E that needs more work to full understood slightly changes 
caused by the obliteration of 00. We are in the way. 

This work, when finished, will imply a sequential work in the kneading 
theory and further study of the entropy of rational real maps. 
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